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Abstract 

' We consider dynamics of the FRW Universe with a scalar field. Using 

Q" 1 Maupertuis principle we find a curvature of geodesies flow and show that 

zones of positive curvature exist for all considered types of scalar field 
{3JTj| potential. Usually, phase space of systems with the positive curvature 

contains islands of regular motion. We find these islands numerically for 
shallow scalar field potentials. It is shown also that beyond the physical 
\ domain the islands of regularity exist for quadratic potentials as well. 

2 ' 1 Introduction 

JO ' It is known from 70-th of the last century that contraction phase of a closed 

C " . Universe filled with a massive scalar field can be followed by expansion for some 



particular initial condition set [TJ [5] • On the other hand, every expansion stage 



. of such Universe is ultimately followed by a contraction one. These two features 

C3 ' of dynamics (which are specific for a closed Universe in contrast to open or flat 

worlds) result in a rather complicated behavior which in some situations can be 
chaotic. 

A chaotic dynamics in massive scalar field cosmology was first found by D. 
Page in [3J, and have been studied in detail in [?J [5] with the corresponding 
discussion on the meaning of chaos in General Relativity. The chaotic dynam- 
ics in question represents an example of a transient chaos with a structure of 
"chaotic repellor" - a countable set of unstable periodic and uncountable set 
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of unstable aperiodic trajectories escaping cosmological singularity. Both sets 
have zero measure in initial condition space. An useful toy model of such kind 
of dynamics (elastic scattering on three discs on a plane) have been described 
in @. 

All these early results have been found for a massive scalar field - scalar field 
with the potential in the form V(</>) = m 2 cf> 2 /2, where m is the scalar field mass. 
Studies of other forms of the potential reveal several different form of dynamics 
with transitions from one form to another (for a short review see [7]). 

The equations of motion can be derived from the General Relativity action 
with a scalar field 

S = J d^xV^^R+lg^d^d^-V^y (1) 

where mp is the constant parameter called the Planck mass, R is the scalar 
curvature of a space-time. For a closed Friedman model with the metric 

ds 2 = dt 2 -a 2 (t)d 2 n i3) , (2) 

where a(t) is a cosmological scale factor, d 2 Q^ is the metric of a unit 3-sphcre 
and with homogeneous scalar field <fi the action (1) gives the following equations: 



a - 



1 \ acf) 2 aV{4>) 



16tt V 2a 2a J 8 4 



- 0, (3) 



4>+^ + V'(<P) = 0. (4) 

with two variables - a scale factor a and a scalar field <f>. 
This system has one first integral of motion 

_ JLm P a(a 2 + 1) + y (f + 2^)) = C. (5) 

This integral plays the role of energy and is equal to zero for cosmological 
solutions of the system (3)-(4). 
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If V(<p) = 0, the system is integrable (in this case a scalar field has the 
equation of state of stiff matter), though the dynamics is finite in time - the 
Universe evolves from Big Bang singularity towards a Big Crunch singularity 
through a point of maximal expansion. There is also another reason why we 
can not develop a perturbation analysis for the system (3)-(5): it can be easily 
checked that multiplying the scalar field potential V by a constant with simul- 
taneous rescaling of the scale factor a and time t does not change the system. 
This means that there is no continuous limit when, for example, m — > for the 
massive scalar field. 

An important property of the constraint equation (5) is that d 2 and <fi 2 enter 
in its left-hand side with opposite signs. Rewriting (5) in the form 



3m 2 P a 2 4> 2 _ 3m 2 p 1 



V{4>) (6) 



8tt a 2 2 8tt a 2 
it is easy to see that 

• There are no forbidden regions in the configuration state. 

• The configuration state is divided to two regions - the region where the 
right-hand side of (6) is positive (and possible extrema of the scale factor 
are located), and the region where it is negative (where possible extrema 
of the scalar field are located). 

The boundary between these two regions is the curve 

2 3 m 2 P 
~ SttVW [<) 

Extrema of the scale factor can exist only in the region where 

2 3 ro| 
" 8tt V(tp) ' 

Using (3) it can be shown that the possible points of maximal expansion (d = 0, 
a < 0) are localized inside the region 

1 m 2 P 



a 2 < 



4tt V{<p) 



while the possible points of transition from contraction to expansion (often 
called bounce) lie outside this region being at the same time inside the region 
of scalar field extrema. 

Zero- velocity points (d = <f> — 0) of a trajectory, it they exist, should lie on 
the curve (7). Numerical studies show that trajectories with these points play 
an important role in the described chaotic structure. In particular, all primary 
(i.e. having one bounce per period) periodical trajectories have zero-velocity 
points as the points of bounce (see numerical examples in [4]). 

Numerical integrations show also that there are regions on the curve (7) 
which can not contain points of bounce. If a trajectory, starting from the 
curve (7) is directed inside the zone of possible extrema of the scale factor, it 
rapidly goes through a point of maximal expansion and evolves further towards 
a singularity. The condition for a trajectory to be directed into the opposite 
zone (the zone of possible extrema of the scalar field) can be written as 

4>/a > d4>(a)/da (8) 

where the function 4>(a) in the RHS is the equation of the curve (7). 

The case of equality in (8) corresponds to a trajectory, tangent to the curve 
(7). This situation was first described in [3], and we call such point as a Page 
point. For the system (l)-(3) the equation for the Page point is [5] 

/ 3 m 2 

V{4> page ) = \j -^.V'ttpage) (9) 

For power-law scalar field potentials the condition (8) is satisfied if 4> > 
4>page > & n d the corresponding part of the curve (7) contains zero- velocity bounce 
points of periodical trajectories. For exponential potentials the condition (8) 
can be violated for all points on the curve (7), and the whole chaotic structure 
disappears [8j . The dynamics in this case can be called regular in the sense that 
the structure of chaotic repellor is absent, and there are no trajectories escaping 
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a cosmological singularity. Global structure of trajectories is the same as for 
V = case - any trajectory starts and ends in a singularity. This does not 
mean that there are no local instabilities of trajectories (see the next section), 
however, as the dynamics is finite in time these instabilities are not important 
for the global picture. 

In the next section we study local instabilities of trajectories, in Sec. 3 we 
describe a new regular regime for the system (3)— (5) , which appears when the 
scalar field potential is sufficiently shallow. 

2 Local instability 

The system (3)-(4) can be considered as a Hamiltonian system with 

2tt 9 1 9 3m 2 P o Tr . ,. o 
H = -^* + l&*--l£- +V M a > (10) 

where the canonical variables are q l = (a, </>). 
We rewrite (10) in the form 

H = a ij p iPj + V (11) 

where 

a J = diag{-- j, — 

V = V((j))a 3 - 3m^a/(87r). 

According to Maupertuis principle we can introduce an auxiliary metric 
space with geodesies coinciding with trajectories of initial system. This method 
of geodesies have been used actively since the middle of the last century [S] in 
various physical problems from Yang-Mills fields [TU] to cosmology [11] . Cal- 
culating the Riemann curvature 1Z of this space we get the criterion of local 
instability in the form 1Z < 0. For the system with a Hamiltonian in the form 
(11), the resulting curvature is [T21 ITB"] 
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n = {n - l) y 



(C _ V) _9 2 V _« (n~6)dVdV 



dq l dqi 4 c?g* dqi 

where n is the number of degrees of freedom. 

Taking C = we get the curvature for cosmological solutions of our system: 



K = 



2 (§I^ a _ a 3y)3 



" / 3 ml „ \ , / ,„ 7r 9 \ 3 ml 1 
V 8 7r / \ rrip / 16 i a 

(12) 

where ' denotes partial differentiating with respect to a scalar held 0. 

The curvature diverges at the curve (7), and this is a consequence of the fact 
that the supermetrics a 13 is not positively definite. This property is typical for 
dynamical systems of General Relativity and discussed in }13j . 

Three different cases are plotted in Fig.l. The first plot (Fig. 1(a)) represents 
the situation of a massive scalar field with the potential V — m 2 </> 2 /2. We 
know that a chaotic repellor exists in the corresponding dynamical system. The 
Fig. 1(b) is created for the potential V = cosh(0) — 1, for which it is known that 
the structure of chaotic repellor is absent. 

We can see that in both cases there are zones with 1Z > as well as zones 
with 1Z < 0. It means that studying only local properties of our system, it is 
impossible to understand its global features. The chaotic system does not belong 
to the class of Anosov's U-systems (this would require 1Z < in all points of the 
configuration space), which posses maximally strong statistical properties [14] . 

However, some qualitative difference in configuration of 1Z < and 1Z > 
zones exists between the chaotic case in Fig. 1(a) and the regular case in Fig. 
1(b). It can be easily checked that the curve 1Z = intersects the boundary curve 
(7) at the Page point (9), so the picture for the system with Page points (as for 
power-law potentials) and without Page points (for steep potentials) should be 
different. Moreover, we can see from the plots that the area of possible bounces 
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(which is located between two hiperbolae) belongs to a stable zone in a regular 
case, and some part of this area is located in unstable zone for the chaotic case. 
This can be checked further for the potential V = (exp(4>/(f)o) 2 +exp(— 4>/c/)q) 2 ) — 

2 with cf>o < 0.96mp, where chaotic repellor exists for intermediate values of 4> 
(Fig. 1(c)) [S]. All this lead to a suggestion that negative curvature in zone of 
bounces is required for the described type of chaotic behavior, though more 
effort is needed for better understanding this connection. 

3 Islands of regularity 

The main goal of this section is to claim that another type (different from 
described in Introduction) of regular dynamics exists for the system (3)-(4). 
This type of dynamics is known for a variety of dynamical systems [15-18]. We 
have found it in the cosmological dynamics (C = 0) with shallow potentials. 
In a general situation this behavior exists even for quadratic potential if the 
constant C from the constraint equation (5) is sufficiently large and negative Q. 
In all our numerical studies we use the units mp/yl&n = 1. 

It has been already reported that in these two cases the picture of chaotic dy- 
namics changes significantly in comparison with the cosmological case of power- 
law potentials. In Fig. 2 we reproduce the plot from [7J, where zones in ini- 
tial condition space leading to at least one bounce is shown for the potential 
V = m 2 <p 2 /2 and negative C. A trajectory starts at the point of maximal ex- 
pansion, so d = 0, and using the constraint equation (5) for eliminating 4> we 
get 2-dimensional space of initial conditions (a, </>). In Fig. 4 the tj> = section 
of this space is shown depending on C. When some bands, being separate for 
small \C\ merge with \C\ increasing, structure of trajectories changes. Trajec- 
tories starting from the wide zone formed by merged bands show complicated 

1 In nonstandard cosmology this type of dynamics exists also in braneworlds [19] and in 
theories with quadratic curvature corrections to Einstein gravity [20]. 
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behavior, which can not be described in a simple way, like in the C — case. 
That is why this type of behavior have been called as "strong" or "less ordered" 
chaos. However, the dynamical behavior in this case appeared to be even more 
complicated in the sense that islands of regular dynamics can be distinguished 
in a chaotic " sea" . 

To see this we have studied an appropriate Poincare section, taking for this 
purpose a point of maximal expansion, i.e. the point with a = 0, a < 0. In 
Fig. 3 the points of this Poincare section are plotted for quadratic scalar field 
potential with m = 1.0 and C — —300. Initial values of the scale factor in 
the point of maximal expansion vary from 7 to 12. Island of regular dynamics 
is easily distinguishable. It should be pointed out that this regular regime is 
of completely different nature than regular dynamics for a steep scalar field 
potential. Now we can see a quasi-periodic behavior, and trajectories from this 
part of the phase space never fall into a singularity. A particular strongly stable 
periodical trajectory exists for initial a ~ 8 (in general, this value is a function 



In the vicinity of this trajectory the points of Poincare sections form rather 
regular ovals, those from more distant orbits show some structure in their 
Poincare set. In the transition zone between regular and chaotic zones non- 
trivial structures can appear, like shown in Fig. 4. Here we can see a Poincare 
set consisting of 9 disconnected ovals. Structure of zone between chaos and 
regularity requires further investigations. 

We have seen similar islands of regularity in the cosmological solutions (i. e. 
with C = 0) of the system (3)-(4) with a family of shallow Damour-Mukhanov 
potentials [3T] studied with respect to their chaotic properties in [22- This 
family has the form 



of C). 




(13) 
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with three parameters - Mq, q and ipo and interpolates between a quadratic 
form of the potential for small <f> and a power-law form V{4>) ~ <$> q for large <^>. 

It was found numerically that the described event of bounce band merging 
exists for a subfamily of (13) with q < 1.24 and tfio < tfio, where <j>o is some 
function of q j22j . In the table below we compare <fio(q) with the value of 4> 
for which the island of regularity appears (it exists for (f>o < 4>q)- They are 
close enough, and it indicates that bounce zone merging is a good indicator 
of regularity in our system. The difference between these two values grows 
significantly only in the vicinity of the critical value of q ~ 1.24 when this type 
of dynamics disappears for any </>q. 
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To our knowledge, this is the first reported case of stable periodical cosmo- 
logical regime in General Relativity (other known stable examples like described 
in [23) require modifications of Einstein gravity, some other solutions periodic 
with respect to scale factor evolution in time are accompanied by a monotonic 
growth of a scalar field [24 ). This regular island is surrounded by a "sea" of 
chaotic trajectories. In general all of them are finite and end in a cosmological 
singularity, though existence of chaotic infinite trajectories of non-zero measure 
is not currently ruled out in our numerical studies. 

4 Conclusions 

We have studied local instability of trajectories representing evolution of a closed 
isotropic Universe with a scalar field using the Maupertuis principle. Zones of 
positive and negative curvature 1Z of the configuration space are constructed. It 
is shown that the system in question is not a U-system, and there are zones with 
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positive and negative 1Z. It is interesting that though zones with different sign 
of 7Z exist for both power-law and exponentially steep potential, in the latter 
case (when the structure of a transient chaos is absent) the configuration of 
zones differs qualitatively from the case of less steep potential (when the chaos 
is present). Our empirical data suggests that a crucial feature for chaos to exist 
is local instability of trajectories in the narrow zones when bounce can occur. 

A typical property of dynamical system with nonnegative 1Z is that their 
phase space contains positive-measured islands of regular dynamics [16] . We 
show that the dynamics of this type exists in the case of sufficiently shallow 
potential as well as in another similar systems with merged bands of initial 
conditions leading to bounces. Trajectories in these islands never reach the 
cosmological singularity, as well as inflationary regime, and have trapped for 
infinite time intervals in a certain zone of small scale factors. 
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Figure 1: Zones of negative curvature (dashed) for the quadratic scalar field 
potential (a), for the potential V — cp^h(0) — 1 (b) and V = (exp((f> / '4>o) 2 + 
exp{—(j)/ fa) 2 ) — 2 with <po < 0.9Qmp (c). Possible bounces can be located 
between bold and dashed lines. 
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Figure 2: Zones of initial scale factor for trajectories which have bounce if 
started from their points of maximal expansion. For large enough negative C 
separate intervals of initial scale factors merge. 
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Figure 3: A Poincare section for trajectories starting in the zone of bounce 
intervals merging. The scalar field potential V{<p) — to 2 </> 2 /2, m = 1.0, the 
constant C = —300. Regular structure around initial a ~ 8, <f> — is clearly 
seen. 
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Figure 4: A Poincare section for a particular trajectory near the boundary of 
the island of regularity. The potential and C are the same as in Fig. 3. 
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